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ABSTRACT 


Calculation  procedures  are  developed  using  axisymmetric 
and  two-dimensional  turbulent  jet  theory  for  predicting  the 
static  lift  performance  of  a two-dimensional  power-augmented- 
ram  wing.  The  resulting  lift  prediction  is  found  to  be  in 
good  agreement  with  experimental  data. 
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INTRODUCTION 

The  static  lift  and  drag  performance  problem  for  a two-dimensional 

power-augmented-ram  wing  has  been  solved  using  incompressible  potential 
1 2 

flow  theory.  ’ However,  to  obtain  good  lift  correlation  with  three- 

4t 

dimensional  experimental  data,  it  has  been  necessary  to  extend  the 
potential  flow  theory  by  including  the  effect  of  turbulent  jet  mixing. 
The  turbulent  mixing  theory  is  used  to  calculate  an  equivalent  potential 
flow  under  the  wing  so  that  the  formulas  developed  from  potential  theory 
may  be  used  to  calculate  lift  and  drag. 


DEVELOPMENT 

In  Figure  1,  the  two  dimensional  jet  is  shown  reacting  with  the 

ground  at  ambient  pressure  P to  provide  momentum  M.  flowing  toward  the 

a 1 


^Gallington,  R.U. , "Sudden  Deceleration  of  a Free  Jet  at  the  En- 
trance of  a Channel,"  DTNSRDC  Departmental  Report  ASED  350  (Jan  1976). 

2 

Galllngton,  R.  W.  and  H.  R.  ChapUn,  "Theory  of  Power  Augmented 
Ram  Lift  at  Zero  Forward  Speed,"  DTNSRDC  Departmental  Report  ASED  365 
(Feb  1976). 

^Reported  informally  by  F.  H.  Krause  and  R.  W.  Galllngton 
("Static  Performance  of  a Power  Augmented  Ram  Wing,"  DTNSRDC  ASED 
TM- 16-76-76,  May  1976). 


ram  wing  with  height  h,  chord  C,  and  flap  height  t2.  Conserving  hori- 
zontal momentum  in  the  ground  reaction 


M = M cos9  + M, 
1 o 4 


(1) 


But  if  total  scalar  momentum  is  conserved  in  the  reaction,  M,  = M - M, . 

4 o 1 

Thus  from  equation  (1) 


= (1  + cos0)/2 


(2) 


With  a turbulent  jet,  it  is  assumed  that  the  velocity  profile  is  "frozen" 
during  the  ground  reaction  and  this  same  result  is  obtained. 


M. 


When  the  duct  flow  under  the  wing  is  filled,  part  of  the  flow  will 
go  over  the  wing  at  angle  0.  Assuming  the  flow  under  the  wing  is  thoroughly 
mixed  so  that  the  pressure  and  y^oclty  are  uniform,  one  may  write  a hori- 


zontal momentum  balance 


P h * (P.-P  )h  = M,  - M.cosB  - pU‘h 
c 5 a 13  5 


(3) 


2 


r 


Assuming  is  uniform,  one  may  write  a horizontal  momentum  balance  under 
the  wing  and  along  the  flap 


P„h  - ujt,  - ujh  + D 
c £.2  3 r 


(4) 


where  Is  drag  due  to  the  flap,  (P-P^)dy. 


Further,  by  assuming  the  flow  under  the  wing  Is  one  dimensional,  one  may 
write  the  energy  equation 


P^h  - JspU2h[l  - (U5/U2)^] 


(5) 


Using  (Uj/U2)  ■ (t2/h)  from  continuity  and  M2  “ pU2t2,  Equation  (5)  becomes 


P^h  * SsM2th/t2  - t2/h] 


(6) 


From  Equation  (4) 


Dj.  “ *sM2[h/t2  + t2/h  - 2] 


(7) 


The  lift  of  the  wing  is  L » Using  Equations  (6)  and  (7),  one  may 


form  simplified  dimensionless  ratios  for  comparison  with  experimental 
data 


Mo 


Mj^C  " ** 


(8) 


D M 

m7  ' M^  [h/tz  + t2/h  - 2] 


(9) 


From  Equations  (8)  and  (9)  It  is  clear  that  the  flap  position  t2  and  the 


ratio  t2/h  substantially  control  the  lift  and  drag  of  the  power-augmented- 


ram  wing  at  zero  airspeed.  As  long  as  t2/h  Is  small,  the  duct  under  the 


i 


3 


wing  remains  filled  and  some  of  the  flow  will  be  lost  over  the  top  of  the 
wing  as  H^cosB,  shown  In  Equation  (3).  As  (t^/h)  1 the  duct  may  become 

underfilled  and  rapidly  lose  lift.  This  occurs  when  the  jet  Is  thin  and 
has  not  entrained  sufficient  mass  to  stagnate  near  the  leading  edge  and 
fill  the  duct.  In  effect,  when  the  duct  Is  underfilled,  the  jet  streams 
under  the  wing  at  ambient  pressure  and  practically  no  lift  Is  generated. 
The  effect  of  turbulent  mixing  Is  Included  by  using  a two  dimensional 
turbulent  jet  velocity  profile  to  describe  the  flow  from  the  jet  origin 
to  the  wing.  Using  ( )'  coordinates  along  and  normal  to  the  centerline, 

U’  = [3t^U^0/4x']’^[l  - tanh^n]  (10) 


n = ay'/x'  (11) 

where  a is  an  empirical  constant,  a = 7.67,  and  t is  the  thickness  of 

o 

the  jet  at  its  origin  where  it  has  uniform  velocity  U'^  = U^.  The  normal 
component  of  velocity  is  neglected  in  these  calculations.  Since  the 
stagnation  streamline  in  the  ambient  pressure  ground  reaction  will  divide 
the  flow  into  and  M^,  one  may  use  Equation  (10)  to  write 

R, /M  =(1/M  )[  pU'^dy'  = (3/4)  [2/3-tanhTi  +(l/3)tanh^n  ] (12) 

T.  o o Jn^  s s 

Substituting  Equation  (12)  Into  (2) 

3 

3tanhn  - tanh  n + 2cos0  = 0 (13) 

s s 

This  equation  may  be  easily  solved  numerically  for  -1 ^tanho^  <.0,  where 

O^B^-ir/l,  by  using  the  Newton-Raphson  method  with  3tanhn  • -2cos0 

®1 

for  the  first  approximation.  Having  solved  for  tanhn^,  one  may  now 

^Schllchtlng,  H. , "Boundary  Layer  Theory,"  4th  ed.,  1960,  McGraw- 
Hill,  New  York,  pp.  605-611. 
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! 
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write  the  equation  for  mass  flux  flowing  toward  the  wing. 


Q,/Q  = (1/Q  )f  pU'dy'  * [3x'/4t  o]"*[l-tanhn  ] (1^) 

1 o o Jn^  o s 


Similarly,  one  finds  the  total  mass  flux  approaching  the  ground,  i.e.. 
Immediately  before  the  ground  reaction  so  that  tanhn^  = -1»  is  given  by 


Q/Q  = [3x'/t  0]*^  » 0.6254[xVt  ]'^ 
o o o 


To  obtain  good  correlation  with  the  experimental  data.  It  was  neces- 
sary to  assume  that  when  t^/h  0 the  flow  under  the  wing  comes  from  the 
flow  with  the  highest  energy,  i.e.,  about  n ^ 0 since  from  Equation 
(10)  this  is  the  maximum  velocity  available.  When  one  considers  that 
the  actual  flow  toward  the  wing  was  generated  by  one,  two  or  three  fans 

during  the  experiment,  where  the  flow  expands  to  form  a quasi-two-dimen- 
4 

sional  jet.  then  it  appears  correct  that  the  greater  inertia  of  the 
flow  about  each  core  should  cause  the  highest  energy  stream  tubes  to 
spread  and  form  a sheet  along  the  ground  with  lower  energy  flow  displaced 
above  in  sheets  of  decreasing  energy.  Thus  the  momentum  and  mass  flux  of 
the  quasi-two-dimensional  jet,  which  is  assumed  for  the  ambient  pressure 
momentum  balance  and  for  calculating  the  flow  under  the  wing,  should  be 
equal  to  the  momentum  and  mass  flux  calculated  at  the  ground  plane  for 
the  axi-s5nnmetrlc  turbulent  jets  from  the  fans. 

The  equations  for  the  turbulent  axisymmetric  jet  calculation  are 


"’a  -(#)■“  »aV«‘  + 


h*  o.r ' /x' . 

A A A 


^Knystautas,  R. , "The  Turbulent  Jet  from  a Series  of  Holes  in 
Line,"  The  Aeronautical  Quarterly,  Vol.  XV  (Feb  1964)  pp.  1-28. 


where  0 is  the  diameter  of  the  Jet  at  its  origin  where  it  has  uniform 
velocity  U'^  “ and  o^  = 15.174.  Again,  the  radial  component  of 
velocity  is  neglected.  Forming  the  axlsymmetrlc  mass  flux  Integral 

^A^^^o  “ 2iTpU^r’dr’  = 0.4566xj^/D  (17) 


Since  the  pressure  is  ambient  throughout  the  jet,  = 

(iv/4)pD^U^  . From  Q/Q  = 0.6254[x'/t  ]**  and  Equation  (17),  the  mass 


flux  ratio  is  taken  as  1 until  x!  = 2.19D  or  x'  = 2.56t  . Beyond  these 

A O 

critical  values,  the  entrainment  process  Increases  mass  flux  until  re- 
stricted by  the  ground.  This  can  be  seen  in  Figure  (1),  where  the  two- 
dimensional  flow,  M^,  can  only  entrain  from  above. 

Using  average  values,  ( ), 


(t/tg)  - (Q/Qq)  /(M/M^) 


is  the  thickness  of  the  jet  since  (Q/Q^^)  = (M/M^)  = 

(U/U^)^ (t/t^) . These  average  values,  when  properly  formed,  may  be 
used  to  calculate  the  lift  and  drag  ratios  given  by  Equations  (8)  and 
(9),  respectively.  For  a two-dimensional  jet  reaching  the  ground 
immediately  ahead  of  the  wing,  the  thickness  will  be 


t = 3x'/a  * 0.391x'  (1 

Following  Equations  (13)  and  (11)  , the  mass  flux  and  momentum  going 
over  the  wing  will  be 


Q3/Q0  * [3x'/4t^a]^(l-tanhii2) 


M,/M  = [2  - 3tanhn-  + tanh  n ]/4 

Jo  J J 


Using  this,  the  mass  flux  and  momentum  flowing  under  the  wing  may  be 
readily  found  by  subtracting 


6 


[3x'/4t  a]  ^tanhn,  - tanhn_) 
O j s 


«2/«„  ■ ("2'”o> 


3 3 

[3tanhii_  - tanh  n - 3tanhn  + tanh  n ]/4  (23) 

^ 3 8 S 


Then,  using  Equations  (19),  (22)  and  (23)  into  (18) 


- 2 3 3 

t-/t  = (tanhn,-tanhri  ] /(3tanhn_-tanh  n-i~3tanhn  +tanh  n ] (24) 

^ J S ,5  ^ 8 s 


If  the  ratio  of  h/t  is  known 


(t^/h)  - (t2/t)/(h/t) 


may  be  used  for  calculating  the  lift  and  drag  ratios  of  Equations  (8) 
and  (9).  As  mentioned  earlier,  the  flow  for  small  values  of  (t2/h) 
should  be  calculated  about  n 0;  thus.  Equation  (24)  Is  used  for 
tanhn,  i |tanhn  | . Up  to  this  condition,  use 

j ^ 8 


tQ2/Qol^  = [3x'/t^altanh^n2 


M,/M  - [3tanhn,  - tanh  n-,1/2 


where  0 < tanhn~  < tanhn  , to  calculate  t„/t. 

z s z 
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The  correct  estimate  of  h/t,  used  In  Equation  (25),  Is  essential 

for  good  correlation  with  the  experimental  lift  data.  Since  the  lift 

Is  calculated  as  If  the  ambient  pressure  ground  reaction  occurred 

Immediately  ahead  of  the  wing  leading  edge,  one  must  determine  an 

equivalent  total  mass  flux  and  resulting  t for  a two-dimensional  Jet 

which  matches  the  conditions  of  the  experiment.  From  Equations  (17) 

and  (18),  and  using  M/M^  » 1 Immediately  before  the  ground  reaction, 

one  may  write  the  equivalent  t for  axlsymmetrlc  jets  with  spacing  b 
2 

assuming  t^  » ttD  /4b, 

t » Tr(0.4566x)^/(4b)  - 0.164x^/b  (28) 

Similarly,  for  the  two-dimensional  jet  one  finds 

t = (0.6254)^x  = 0.391X  (29) 

These  two  equations  for  free  air  mixing  are  plotted  In  Figure  (2)  to 

give  a simplified  model  of  the  entrainment  process.  The  value  of 
2 

irD  /4b  for  t^  corresponds  to  Q/Q^  “ 1 and  Is  taken  constant  until 
X 2.19D  where  the  axlsymmetrlc  entrainment  formula.  Equation  (28), 
starts  to  Increase  t.  From  this  point,  t Is  given  by  Equation  (28) 
until  X > 2.38b  where  the  axlsymmetrlc  jets  have  spread  Into  a quasl- 
two-dlmenslonal  jet  and  equation  (29)  must  be  used.  If  the  jet  reacts 
with  the  ground  ahead  of  the  wing,  some  of  the  free  air  entrainment 
will  be  Inhibited  during  the  ground  reaction  so  that  t^  will  only  con- 
tribute about  0.8  Its  calculated  value.  Similarly,  the  entrainment 
during  the  ground  run  must  be  taken  at  0.5  Its  calculated  value  since 
the  two-dimensional  jet  along  the  ground  can  only  entrain  from  above. 

With  the  estimate  of  t from  Figure  (2)  and  h from  the  geometry 
of  the  experiment,  one  may  calculate  t2/h  and  various  values 

of  n using  Equations  (18)  and  (22)  through  (27)  and  In  turn  calculate 
a plot  of  Lh/Mj^C  versus  t^/h  from  Equation  (8) . Note  that  equations 


Figure  2 - Suonary  of  Turbulent  Mixing 


(2)  and  (13)  must  also  be  used.  An  example  Is  shovm  In  Figure  (3). 

A direct  calculation  procedure  may  be  devised  by  assuming  a two 

dimensional  turbulent  jet  velocity  profile  at  the  wing  leading  edge, 

with  u at  the  surface  or  centerline.  Then 
max 

2 2 

U ■ 1 - tanh  n,  “ sech  n,  (30) 

U 

max 

approximately  describes  this  velocity  profile  ahead  of  the  wing 
where,  using  Equations  (11)  and  (19), 


n 


1 


Using  Equation  (30),  the  momentum  Is 


M. 


PU^  ^1  [tanhn,  - 1 tanh^n.Jn 
max  1 10 


(31) 


(32) 


This  must  match  the  momentum  given  by  Integrating  from  n • 

X s 


2 - 1 3 * 

H,  • pU  t [tanhn  - tanh  n] 

X indx  ^ ^ n 

3 s 


Using  Equation  (13),  one  finds 


tj^  “ (1  + cos  G ) t 


(33) 


(34) 


Following  Equation  (32),  the  momentum  under  the  wing, 

(^)] 


will  be 


(35) 


Substituting  these  results  In  Equation  (8) 


10 


should  be  used  to  avoid  the  nonmixing  region  shown  In  Figure  (2b) 
and  In  turn  avoid  calculating  a dynamic  pressure  at  the  wing  greater 
than  Qq.  These  results  are  summarized  In  Figures  (3)  and  (4). 

Note  that  the  drag  plot  In  Figure  (3)  only  uses  one  half  the  theoretical 
drag  prediction.  This  Is  attributed  to  thrust  recovery  at  the  wing 
leading  edge  which  Is  not  considered  In  the  flow  model. 


CONCLUSIONS 


The  static  lift  performance  for  a power- augmented- ram  wing  can 
be  accurately  predicted  by  using  two-dimensional  and  axlsymmetrlc 
turbulent  jet  theory  to  extend  the  existing  potential  flow  analysis. 

The  drag  prediction  using  this  theory  Is  quite  conservative  and  requires 
further  Investigation  to  obtain  good  correlation  with  experimental  data. 


11 
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Figure  4 - Comparison  Between  Theoretlcal/Experlmental  Results 


